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A UNIFIED FORM OF LAMBERT’S THEOREM 


by 

E. R. Lancaster and R. C. Blanchard 
Goddard Space Flight Center 


INTRODUCTION 

Lambert’s problem, as it arises in most applications, is concerned with the determination of 
an orbit from two position vectors and the time of flight. It has important applications in the areas 
of rendezvous, targeting, and preliminary orbit determination. In this paper a unified form of 
Lambert’s theorem will be presented which is valid for elliptic, hyperbolic, and parabolic orbits. 

The key idea involves the selection of an independent variable x and a parameter q such that 
the normalized time of flight T is a single- valued function of x for each value of q. The parameter 
q depends only upon known quantities. The problem then is to find x for given values of q and the 
time of flight. For less than one revolution, T is a monotonic function of x for each value of q. Thus 
it is an easy task to design an algorithm for finding x. For multirevolution cases, T(x) has a single 
minimum for each q. 

This idea was presented in a previous paper (Reference 1), where a unified formula was given 
for the computation of T from x andq. Detailed sketches were given for T(x, q) , and a simple 
formula developed for the time derivative of the magnitude of the radius vector at the initial time 
in terms of x and given quantities. 

The present paper expands upon the previous one by: 

1. giving complete derivations which were only sketched before; 

2. giving a careful derivation of the classical form of Lambert’s equations, including the 
multirevolution cases; 

3. deriving a number of useful auxiliary formulas, e.g., for the semilatus rectum and for the 
velocity vectors at the initial and final times. 

THE PROBLEM 

Suppose a particle in a gravitational inverse- square central force field has distances r 1 and 
r 2 from the center of attraction at times t x and t 2 . Let c be the distance and 9 the central angle 
between the positions of the particle at the two times, where 0 < 9 < 2tt. 
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Lambert's problem is that of finding the semimajor axis or some related quantity for the orbit 
of the particle, given t 19 r 19 t 2 , r 2 , and#. When Lambert's problem has been solved, other quanti- 
ties associated with the orbit are easily found, as will be later discussed. Using the law of cosines, 
we can express c in terms of t 19 t 29 and Q\ 


C 2 - Fj 2 + r 2 2 - 2fj r 2 cos 0 . (1) 

We define 

G = universal gravitational constant, 

M - mass of attracting body, 

M = GM, 

a = semimajor axis of transfer orbit, 
e = eccentricity of transfer orbit. 

We will follow the common sign convention for a, i.e., a > 0 for elliptic orbits and a < 0 for 
hyperbolic orbits* 

Definitions of other symbols will be given as they are introduced. 


THE CLASSICAL FORM OF LAMBERT'S EQUATIONS 

The path of a particle in an inverse-square central force field is an ellipse, parabola, or 
hyperbola. With origin at the center of attraction, we have, for elliptic motion, 

r i = a(l -ecostfj) , ( 2 ) 

r 2 = a (l ~ e cos <t > 2 ) > (3) 

n(t 1 -i p ) = " e sin <j > 1 , ( 4 ) 

n ( t 2 ~ t P ) = ~ e sin <t 2 , ( 5 ) 

where <P 1 and 4> 2 are the eccentric anomalies at times t % and t 2 , t p is the time at pericenter, and 

n = ( M /a 3 ) 1/2 . 
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If €j is a unit vector pointing towards periapsis and e 2 is a unit vector in the plane of motion 
90° ahead of 7, in the direction of motion, then for the position vectors 7, and 7 2 at timts 1 1 and t 2 

fj = a^cos^j-ejej + a (l - e 2 ) 1/2 (sin <p t ) e 2 , 
r 2 = a(cos <£ 2 - e ) e , + a(l - e 2 ) 1/2 (sin <j > 2 ) e 2 . 

Substituting these equations in 

c 2 = r, 2 + r 2 2 - 27 j • r 2 , 
where the dot indicates scalar product, we have 

c 2 = a 2 (cos 0 2 -cos^ij) 2 + a 2 (l - e 2 )(sin<£ 2 -sint£,) 2 


- 4a 2 [l - e 2 cos 2 -J (<£j + 0 2 sin 2 ^ (^2 “ ^ 1 ) ' 


( 6 ) 


Adding Equations 2 and 3 , 


2a [l - e cos -j (<P l + <t> 2 ) cos -5 (<£ s - <P t )] 


(7) 


Subtracting Equation 4 from Equation 5, 


n(t 2 - tj ) - 4> 2 - <f> 1 - 2e cos 7 + 0 2 ) sin 2 (^2 ) 


( 8 ) 


Equations 6, 7, and 8 determine the three unknowns a, <t> 2 -<t> l , and e cos ( 1/2) + 4> 2 ). Let 


cos 2 (a- + /3) - e cos ~2 [$1 + 4> 2 ) < 0 ^ a + ,8 < 2tt , 


a - fi - <f> 2 - 4> l ~ 2nt77 , 0 a - /3 < 277 , 


(9) 

( 10 ) 


where m is the number of complete circuits made by the particle between times and t 2 . 
Equations 6, 7, and 8 become 


c/2a - sin‘2'(a+/3)sin2( a_ /^) • 


(ID 
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(12) 


( r i + r 2 )/ 2a 


1 - cos ~2 (a + fi) cos (a ~ fi) , 


n (t 2 “ t x ) - 2rr rrr + a - fi - 2 cos (a. + /3) sin ( a " fi) 


(13) 


(0, 2u) 


(0, -2tt 



(0, 0) 


The two inequalities for Equations 9 and 10 are 
geometrically equivalent to the shaded region of 
Figure 1, from which it is evident that 0 ^ a < 2 77 
and -7T < fi < 77 . We can also obtain 0 ^ a < by adding 
the inequalities for Equations 9 and 10; and, if we 
add fi - a to each part of the inequality for Equation 9 
and divide the result by 2, we obtain - (a - fi)/2 < fi<rr 
+ (a - fi)/ 2 , Or - 77-/3 < 77 . 

With appropriate trigonometric identities, Equa- 
tions 11 and 12 become 

cos fi - cos a - c/a , 

cos fi + cos a r 2 “ (r x + r 2 |/a • 

Solving these two equations, 


Figure 1 . 
where we have defined 

S - ( r i + r 2 +C )/ 2 • 

E = s/ 2a , 

K - 1 - c/s . 

Since cos a = l - 2 sin 2 (a/2) , Equation 14 can be changed to 

E = sin 2 (a/2) , 0 — a < 277 . 

Similarly, Equation 15 becomes 

KE = - sin 2 (fi/ 2) , - tt < fi < rr , 

K = (s-c)/s = ( r x + r 2 - c )y/2s 


cos a - 1 - s/a - 1 + 2E , 

cos/3 = 1 + 2KE , 


(14) 

(15) 


(16) 


(17) 


[( r i +r 2) 2_c2 ]/ 4s2 ' 
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Introducing Equation 1, we have 


K - ( T 1 r 2 /2 s 2 ) (1 + cos 0) 


“ ( r 1 r 2 A 2 ) COs2 (^/ 2 ) • 


Substituting Equation 16 in Equation 17, 

sin (yd/2) = q sin (a/2) , 


where 


- 7T — (3 < Tf , 


q - ± VK - ^(rj r 2 ) 1/2 / s ] cos (0/2) . 

Note that the sign of q is taken care of by the angle 0 : 

1 > q > 0 if 0 < 0 < tt , 

0 ^- q > - 1 if t [ < 0 < 2 tt . 

We can introduce E into Equation 13 , since 

n(t 2 -ti) - ( w/a 3 ) 1 2 ( t 2 “ 1 1 ) = (-E) 3 ' 2 T 

where 


T - (8/Vs) 1/2 (t 2 - tj )/s . 


T 


( -£)“ 3/2 




This can also be written as 

T = ( - E) 3/ 2 [2m7T + a- /3-(sina-sin/3)] . 

Substituting Equation 16 into Equation 22, 

T sin 3 ( a/ 2 ) = 2mrr + a - /3 ~ sin a + sin/3 . 


(18) 


(19) 


(20) 

( 21 ) 


( 22 ) 


(23) 
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Equations lfi and 23 with 0 1 a < 2w are Lambert's equations for elliptic motion. Given T and 
q, they are to be solved for a and p t after which it is a simple matter to find all other quantities 
associated with the orbit. 

It is customary in the literature (e.g., Reference 2) to consider T as a function of E (or a) and 
break the elliptic case of Lambert's theorem into four cases, depending upon the sign of q and 
Whether a/2 is taken in the first or second quadrant in Equation 16. The choice of E as the inde- 
pendent variable makes T a double-valued function. This problem can be avoided by choosing a as 
the independent variable. However, an even better choice will lie discussed in the next section. 

By a derivation very similar to that for the elliptic case, one finds for the hyperbolic case 


T E' 3/2 [y - S - ( sinh y- sinh S)j 


(24) 


E = sinh 2 (y/2) , 


(25) 


sinh ( S/2) - q sinh (y/2) 


(26) 


When m = 0, Equations 22 and 24 break down for E = 0 and suffer from a critical loss of sig- 
nificant digits in the neighborhood of E = 0. To remedy this, Equation 22 is written in the form 


T = cr(-E) - qKcr(-KE) 


(27) 


where 


ct(u) - 2 [a res in u 1/2 - u 1/2 (1 - u) 1/2 J/u 3/2 . 

Replacing arcsin u 1/2 and (l -u) 1/2 by series (Reference 3) with 01 a. < n f 

oo 

cr(u) - 4/3 + y ’ a n U n , | U | < 1 , 

n = 1 

a n = 1 * 3 * 5 • • • ( 2n ~ 1 )/ 2 n ” 2 ( 2n + 3 )n ! 

A similar procedure produces the same series for the hyperbolic case. For the parabolic case 
we have E = 0, in which case the series gives 

T = (4/3) (l -q 3 ) . 

Thus with m = 0 we have a series which is valid for elliptic, hyperbolic, and parabolic transfer 
provided |E| < 1 (with 01 a < 77 for the elliptic case). 
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A UNIFIED FORM OF LAMBERT’S EQUATIONS 


As mentioned in the previous section, T is a single- valued function of a. However, a better- 
behaved function is obtained if we choose as the independent variable 

X = cos (a/2) , -1 < X ^ 1 t 

= cosh (y/2) , x > 1 . 

We then have, for elliptic, hyperbolic, and parabolic transfer, 

E = x 2 - 1 ; 

if for parabolic transfer, we let x = 1. 

For the elliptic case, let 

y - sin (a/2) = (~E) 1/2 , 

z = cos (/3/2) = ( 1 + KE) 1/2 , 

1 

f - sin j (a-ft “ y(z-qx) , 

g = cos y ( a “/3) - xz - qE , 

1 

h - ~2 ( s in a - sin /3) - y(x-*qz) , 

\ - arctan(f/g) , 0 < \ < tt . 

It then follows from Equation 21 for the elliptic case that 

T = 2(nT77 +\~h)/y 3 . 

For the hyperbolic case, let 

y = sinh (y/ 2) = E 1/2 , 

z = cosh (y/ 2) - (1+KE) 1/2 , 
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f = sinh-^(y-S) - y(z-qx) , 
g = cosh \ (y - 8) = xz - qE • 

Note that 0 < y - 5 < 00 , since 0 ^ f < 00 . Let 

h = \ ( sinh y - sinh 8) = y(x-qz) . 

It follows that 

-j(y-8) = arctanh ( f/ g) , 

= \ in [( f + g)/(g - f )] . 

= \ In [( f + g) 2 /(g 2 “ f 2 )] 

= In ( f + g) . 

Thus, for the hyperbolic case, 

T = 2[h-Cn(f + g)]/y 3 • 

It is now apparent that, given q and x, the following steps produce T for all cases 

1. K = q 2 

2. E = x 2 - l 

3. If x is near 1, compute T from Equation 27 

4. y = I e| 1/2 

5. Z = ( 1 + KE) 1/2 

6. f = y(z _ q x ) 

7. g r XZ - qE 
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8. If e < 0, \ = arctan (f/g), d = mn + 0 < K <n 

KE>0, d = 4n(f + g) 

9. T = 2( X - qz - d/y)/E. 

The following formula for the derivative holds for all cases except for x = 0 with K = 1, and 
for x = 1: 


dT/dx = ( 4 - 4qKx/z - 3xT)/E . 

If x is near 1, differentiate Equation 27 to obtain 


dT/dx = 2x [qK 2 cr* (~KE) ~cr' (-E)] 


a' (u) - do/du 


2^ na n un_i 


The derivative in the case of x = 0 with k = 1 will be discussed in the next section. 


AUXILIARY FORMULAS 

In this section we will show how to obtain a number of useful quantities associated with the 
two-body orbit, assuming Lambert T s problem has been solved for x. 

In the derivation of Lambert’s equation for the elliptic case, * and / are defined in such a way 

that 

0 2 - - a - 3 + 2xm - 2(\ + nw) . 

From Equation 22, the eccentric anomaly difference can also be written in the form 

<t> 2 “ <t> x ~ (~E) 3/2 T + sin a. - sin 


- y 3 T + 2y( x “ qz) 


(28) 


From Equation 28 we have 




2y( z “ qx) (xz - qE) 


- 2y(x-qz) + 4y 3 q(z-qx) , 
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( 29 ) 


<t> 2 . ~ <t>i ~ sin (<t > 2 - «£j ) " y 3 T - 4y 3 q( z - qx) , 
1 - cos (<£ 2 - <£j ) = 2y 2 ( z ~ qx ) 2 . 


(30) 


We now obtain a formula for the time derivative r at time t t . Kepler's equation in the elliptic 
case can be written in the form (Reference 4) 

W a3 ) 3/2 (t 2 " t i) r + r i ' T i [l -cos (<t> 2 -^ 1 )] / /(Ma) 1/2 - (l - r j/ a ) s in K '^l ) • 

Substituting 1/a = 2y 2 /s , t 2 “ t i = s 3/2 t/(8m) 1/2 , and making use of Equations 29 and 30 gives 
(2/us) 1/2 (z * qx) Tj r t = 2q - (2r iy /s ) (xz - qE) • 

Multiplying through by z + qx, we have, since ( z - qx) ( z + qx) = l - K = c/s and ( z + qx) ( xz - qE) 
~ x + az , 

( 2/{is) 1/2 c r 1 t j = 2qs(z+qx) - 2r 1 (x+qz) 

- 2qz(s-rj) + 2x^Ks-r 1 ) ; 

Ks - r 1 - ( 1 - c/s) s - f j - s - c - Tj = r 2 - s . 


Thus we have finally that 


- (2ms) 


1 / 2 


[qz(s -r,)-x(s- r 2 )]y4r , 


(31) 


At time 1 2 we find 


r 2 - (2ms) 1/2 [x(s " ) - qz(s - r 2 )Jy4 


(32) 


By a similar procedure we can show that Equations 31 and 32 hold also for the hyperbolic and 
parabolic cases. 

Having x, we can find the semimajor axis a, or its reciprocal, 


1/a - 2y 2 /s 
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We know that 


e cos 0 - 1 - r/a , 

esin<£ ~ r f/( fta ) 2 , 

where r is the magnitude of the position vector, and 4> is the eccentric anomaly, at time t. 
Thus for the eccentricity we have 

e 2 = (l - r/a) 2 + (rr) 2 //xa . 

For the semilatus rectum p we have 


p - a (l ~ e 2 ) - 2 r ~ r 2 /a - (rf) 2 //i , 

and for the value of r at the point of closest approach to the center of attraction we have 

= P 

r p 1 + e 


For the speed we have 


V 
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For the component of velocity perpendicular to a radius vector, 


~ Op) 1/2 A - 

The above formulas hold for any type of two-body motion. Either of the subscripts 1 or 2 can 
be placed on r, f , and v. 

If t x < 0 and r 2 > 0, or if r 1 and f 2 have the same sign with 0 > n 9 periaps is passage will occur 
between times t x and t 2 , in which case it may be of interest to compute r p . 

If the plane of motion is known, the velocity v at either time t x or time t 2 can be written in 
any convenient coordinate system, since the components r and v e are known. The plane of motion 
can be found from the position vectors r 1 and r 2 at times t 1 and t 2 , provided they are not parallel. 
Since this is a common case, we will express v 1 and v 2 (the velocities at t x and t 2 ) in terms of r 1 
and r 2 . 



For the velocity V 1 we have 


where v rJ is along and v 8l is in the plane of motion perpendicular to r t and in the direction of 
motion, i.e., in the direction of increasing true anomaly. We have 

? rl = ' 

V $1 = Cj 7, + C 2 r 2 , 

where c 1 and c 2 are to be determined. 

? 1 * V*! = 0 = C 1 r i 2 + C 2 ? 1 ' 7 2 


7 2 • - r 2 v d l sin5 - C 1 • r 2 + c 2 r 2 2 . 



Figure 2. 


Since • r 2 - r x r 2 cos (9 9 we have 

r i c x + ( r 2 cos 6 ) c 2 = 0 , 

(rj cos^jcj + r 2 c 2 ~ v 0 sin 9 . 

Solving for c l and c 2 ,we have 

Vi - (f I -v Jl cot«)(r 1 /r 1 )+ (v fll csc<9)(r 2 /r 2 ) . 
In a similar way, we find 

v 2 = "(v ff2 csc(9)(? 1 /r 1 ) + (f 2 + v e2 cot 0 )(? 2 /r 2 ) . 

Figure 2 shows T as a function of x. Note 
the discontinuity in the slope for x = 0 with 
K = 1. For K = 1, z = | x | . Thus we are led 
to consider four cases: q = ±1 with x > 0, 
q = ±1 with x < 0. Examination of the formulas 
for dT/dx in these cases reveals that if q = 1 
we have a left-hand derivative of -8 and a right- 
hand derivative of 0 at x = 0. If q = -1 we 
have a left-hand derivative of 0 and a right- 
hand derivative of -8 at x = 0. 
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With further analysis we find that the cases (q = 1, x < 0) and (q = -1, x > 0) represent recti- 
linear orbits. 

For m = 0, T is a monotone function of x, making possible a simple numerical procedure for 
solving Lambert’s problem. Figure 2 is for the elliptic case, Figure 3 for the hyperbolic case, 
the parabolic case occurring at x = 1 in both figures. Figure 4 shows a small region of Figure 2 



Figure 3. 


where d 2 T/dx 2 is negative. If the Newton- 
Raphson method is being used to find x, a 
switch should be made in this region to the se- 
cant (regula falsi) method. 

No solutions of Lambert’s problem exist 
in the shaded regions of Figures 2 and 3. 
x = 1 (m > 0) and x = -1 are vertical asymp- 
totes. T - 0 as x * 
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